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Abstract 

We will prove certain general relations in Matrix Product Ansatz for one dimensional 
stochastic systems, which are true in both random and sequential updates. We will 
derive general MPA expressions for the currents and current correlators and find the 
conditions in the MPA formalism, under which the correlators are site-independent or 
completely vanishing. 
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1 Introduction 



One of the fruitful techniques for the study of stochastic systems on one dimensional lattices 
is the Matrix Product Ansatz (MPA)0|0|, which is a generalization of the simple product 
measure, where the steady state probabilities are represented by matrix elements or traces 
of product of appropriate operators. This ansatz when applied to processes with random 
sequential updates, states that the steady state of a process governed by a Hamiltonian of 
the form 



JV-l 



H = h x 



H K,k+1 + h N- 



can be written as: 



IP >= 
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where is a normalization constant. Here A is a column matrix with operator entries 
acting on some auxiliary space F, that is 



A = J2A 



% >- 
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and < W\ and \V > are two vectors in F* and F respectively. Here we have assumed that 
the Hilbert space of states of the chain is 7i = h® N , where h is the p+1 dimensional space 
of one site. The space h is spanned by the vectors \i >; i = 0, • • -p where |0 > stands for a 
vacant site and \i > stands for a site occupied by a particle of type i. Note that we use the 
same symbol | > for a vector in h and F, hoping that no confusion will arise. 
The conditions for stationarity of (El) are the following 0] 



h B A A = X ® A- A® X 

{h N A-X)\V >= 
< W\(hiA+X) = 
where X is a suitably chosen vector with in general operator entries, i.e: 

x 2 



(4) 

(5) 
(6) 



x = Y,x i \i>-- 



i=0 



(7) 
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When applied to processes with backward sequential (BS) updates in discrete time with 
updating operator 

T_ T rpBrpB rpB rp / q\ 

— J 1 J 12 J 23 " " I N-1N I N yo) 
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where T\ and are the boundary terms and the rest of operators implement the bulk 
dynamics, the steady state \P > is written in the same form as in (^j, but the MPA realtions 
are replaced by 

(9) 



T B A <g> A = A <g> A 
T N A\V >= A\V > 
< W\T X A=< W\A 
where A is a new operator-valued vector in h. 

A 2 



(10) 
(11) 



A = j2Ai\i >- 



(12) 



V a p J 

The form of the relations for forward sequential (FS) update is obtained from (|9"|-|iTD by 
interchange of A and A. 

This ansatz has been applied to the Asymmetric Simple Exclusion Process (ASEP) under 
various circumstances and with different kinds of updates ( see || ,|J and and references 
therein.) To explore the usefulness of MPA beyond the above simple cases, there have also 
been a number of attemps to study the algebras associated with more complicated pro- 
cesses || P| fllCl |TlI| and fl2|| , or conversely, the processes associated with more complicated 
algebras 1 113| 1 1|14|| 1 15| ||16|| and |T7|]. In pursuing the latter line, that is begining from algebras 
and searching for processes, one usually notices that the kind of algebra severely restricts 
the kind and the rates of the processes, especially on open systems. 

The aim of this paper is to start from the general MPA algebras (§-§|) and (p|~pTTj) and derive 
some general relations, which are valid for a large class of processes. These relations put 
some general constraints on the algebra that one wants to start with for the formulation of 
a process in the MPA formalism. We first state our main results the proof of which will be 
found in the forthcoming sections. 

Main Results: 

We consider the family of processes on an open chain, described by Hamiltonians of the type 
(]!]), where there are a number of species of particles interacting in the bulk. The processes 
in the bulk are quite arbitrary, i.e. particles can be created, annihilated or can coagulate or 
decoagulate. However we assume that some species of particles are conserved, so that for 
each conserved species say the z-th one, a conserved current J % can be defined. Particles are 
injected to or extracted from both the left and the right ends and a driving force may also 
be present. 

Due to this arbitrariness we do not use the detailed form of any particular Hamiltonian or 
the associated algebra, or any of its representations therof, but only use the basic formulas 
of MPA. Our main results are the following: 

1-If a species say the i-th one is conserved, then its current is given by: 

< wiXi^-^v > 



< J 1 >-- 



< W\C N \V > 



2 



where in the ordered updates Xj = A t — Aj. 

2- If Y^i Xi = , then the n-point correlators of all the conserved currents are independent of 
the distances between the points. This is true even in finite open chains and is independent 
of what the other species are doing. 

3- If a set of {X^X^, • • •} are c-numbers, then the connected correlation functions for the 
currents of particles {ii,i2,---} vanish in the thermodynamic limit. Moreover in a finite 
chain with N sites, one has for the above set of particles: 

< j r j s > N =< rr > N =< j r > N < j s > 7V _i (13) 



< .rrf > N =< J r > N < r > N -i< f > N -2 (u) 

etc. 

Remarks: 

a:The above properties when supplemented with the results from numerical computations 
or simulations for small chains, may give us an idea of what kind of algebra can or can not 
be used for an exact solution of the problem in a large chain. 

b:The above results are particular to open systems and do not apply to processes on closed 
systems. 

4-Under any circumstances the currents and densities of the conserved species in backward 
and forward sequential updates are related as follows: 

< n* >^ (k)— < n { >^ (k) =< >=< JV > 

where the direction of the arrows indicate the type of update. This result is a generalization 



of a previous one first derived in [19 



Remark: In deriving the above results we have not used any particular representation of 
the MPA algebra, i.e: like the Fock space representations of fl9| or ||. 



2 General relations in MPA algebras 

In order to deal with all the different updates in a uniform manner we define the following 
objects 

hf.^l-Tt, h N :=l-T N , h B :=l-T B (15) 
and define in the BS and FS updates the vector 

X:=A-A. (16) 

With these defintions, the MPA relations for these ordered updates take the following 
form, compared with those for the random update (see eqs. (1-4)). 

h B A® (A-X) = X®A-A® X (17) 

(h N A- X)\V >=0 (18) 
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<W\{h 1 A+X) = (19) 



Note the similarity between the equations for different updates. In fact (|5]) and (JL8|) are 
exactly the same. This rewriting could be of no use were it not for the fact that usually (see 
for example [2I| ||18|| ), the evolution operator T B of a process in ordered update is, modulo 
a redefinition of parameters, nothing but 1 + h B where h B is the Hamiltonian of the same 
process in random sequential update. Therefore this rewriting allows one to immediately 
write down the MPA relations for any of these updates, once they are known for one of 
them. Furthermore, in those situations where a map may exist between the MPA algebras 
for different updates, such rewriting facilitates the search for such a map. We will see an 
example of this later on. Before proceeding further we fix some notations and conventions. 
In the Hilbert space h we define a reference state < s\ as 

i=p 

<a|:=X;<i| = (l,l,...l) (20) 

i=0 

The reference state for the space h £g> h is defined as < ss\ :=< s|<8> < s\ and similarly for 
tensor products having more factors. This state is used to write the sum of entries in a 
column of a vector or a matrix in closed form. The following operators acting on the space 
F are also useful: 

i=p i=p 

C :=J2 A i=< s\A> K := Y^ x i =< s\X > (21) 

i=0 i=0 

The normalization constant is given by =< W\C N \V >. Conservation of probabilities 
imply that we should have in all three updates 

<: ss \h B =< s\h =< s\h N = (22) 

Actually in ordered sequential updates, conservation of probability need not hold in every 
single update, but only in one complete update. Therefore one should only have < S\T =< 
S\. In writing ( |2"2"D we are assuming conservation of probability at each single update, i.e. 
< ss\T B =< ss\. The known sequential updating procedures for the ASEP fall within this 
class. 

Multiplying both sides of the MPA relations (|4] and (5|) or (|17| and |1^) (from the left by < ss\ 
or < s\ where appropriate and using (pl|) and (p2|) we obtain the following relations which 
are valid in all kinds of updates 

K\V >=< W\K = [K,C} = 0. (23) 



This relation has been obtained first in [I9~|. 

We now consider a processes in which some of the particle species are conserved locally. In 
this case the local operator h B k+1 does not change the number of particles in a conserved 
species say the i— th one, on the couple of sites k and k + 1. Thus we have 

irt + H + i,hf k+1 } = (24) 
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where f l is the number operator of particles of type i, i.e. f l \j >= 5ij\j >. Applying the 
state < ss\rl + r l k+1 on both sides of (f| or [T7]) and using (^T] and |22|) we obtain: 

[X l ,C} = [A l ,K] \/i (25) 

This relation holds for each conserved species seperately and is a new constraint on the MPA 
algebra for such processes. 

Remark: In a process consisting only of exchange of particles, more detailed relations can 
be obtained. However eq. ( |2"oD is valid independent of the interaction of the other species. 



3 The MPA expression for the currents 

The average density of particles of type i at site k is defined as < r % k > (t) :=< S\fl\P(t) > 
where < S\ is the reference state of the whole lattice. We will now obtain general MPA 
expressions for currents of conserved species of particles. 

At least in the RS update, the MPA expressions for the current of a conserved species can 
be derived simply by calculating the current at one of the boundary links, say the rightmost 
link. This current is the product of density at the rightmost site and the extraction rate 
of that species. However we prefer to follow a different approch and calculate the currents 
directly in the bulk, in order to calculate also the current correlators (see section 3.1). 
1) The RS update: 

According to the Hamiltonian formulation of Markov processes, we have 

j t <Tl>=<S\[H,rl]\P(t)>. (26) 
Due to the form of H this can be rewritten as 

j t < rl >=< S\ [ht hk + h* k+1 , fl]\P(t) > . (27) 
Combination of (^) and ( |2"7D now yields 

j t < ri >=< suht^nwm >-< s\[h» k+1 ,ti +1 ]\p(t) > (28) 

which can be written as a continuity equation: 

j t < rl >=< 4 > - < T k+1 > . (29) 

with the current of particles of type i into site k being 

<4>=<S\[ht ljk ,Ti]\P(t)> (30) 
To find the MPA relation for the steady state currents we rewrite (|30|) as 

<J*>=<W\Ji\V> (31) 
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where J k is the following operator acting on F 

1 



Jl = — < S\ [h B _ 1}k , rl]\A®A®---A>. (32) 

Note that here we are looking at \A® A® ■ ■ ■ A > as an operator valued vector in TC. Using 
(gj), we find 

J l k = -^-C k - 2 < ss\[h B , 1 <g> r l ]\A <g> A > C N ~ k =: -^C k - 2 M l C N ~ k (33) 

The operator Mi is calculated as follows: 

Mi = < sslh 3 ^®^) - (l®r i )h B \A®A> 
= - < ss\l®r i \X®A-A® X > 

= -KA + CXi (34) 

where in the second line we have used (|j and |22D and in the third we have used (|3|) , (0) and 
(H). Thus using (H), we find: 

p < w\c k -^c N - k \v > 
< J >_ < ty|c^|y > (d5j 

We can now use (p5|) and rewrite CXj = XjC + KAi — and move the two K's to the 
left and right where their action on the vectors < W\ and \W > vanish. In this way we can 
move Xi completely to the left and obtain: 

Tl < w\x,c- l \v > 
* J >= < w\c«\v > (36) 

This kind of relation for conserved currents has already been obtained for specific 2-species 
algebras[n] where all the particles are conserved. Here we derive it in quite general form. 
2) The BS and FS updates: 
In the BS update we have: 

< rl > (t + 1)- < rl > (t) =< S\rtT\P(t) >-< S\H\P(t) > (37) 

Using the property < S\T =< S\ , we can write (|37| ) in the form: 

< tI > (t + 1)- < Tl > (t) =< S\[Tl,T]\P(t) > (38) 

Taking the structure of T into account (see (§)), using the notation T k B t := T k B k+l ■ ■ ■ T B l+l ( 
with T B N+1 := T N ) and the property < S\T k B =< S\, we can rewrite eq.(|3~8"D as: 

<rl>(t + l)-<Tl>(t) = <5|[f^Tf_ ljfc T^ +1 ]T4 ljJV |P(*)> 

= < S\{il,T k \ k ]T k B N \P(t) > 
+ < S\[fl,T B k+1 ]T B +hN \P(t) > (39) 

where in the last line we have used < ss\T B _ l k =< ss\. Using the local conservation law 

[rl + fl +1 ,T B k+1 }=0 (40) 
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Equation ( p9|) can be written in the form of a continuity equation 

< rl > (t + 1)- < rl > (t) =< 4 > (t)- < 4 +1 > (t) (41) 

Where J\ is the number of i-particles which, in the interval between t and t + 1, leave site 
k — 1 and enter into site k, and is given by 

< Jl >=< S\[fl,T k \ k }T k B N \P(t) > (42) 

The MPA relation for the current < J l k > is now obtained along the same lines as in the RS 
update, namely 

< 4 >=< W\J l k \V > (43) 

where 

4 =< S\ [rl T k \ k ]T B N \A ®A®---A>. (44) 
Using (Q) and ( pT|) we obtain 

4 = < S\[r k ,T k B _ ltk ]\A®A®---A®A---A> 

k-l,k 

= C k - 2 < ss\[1®t\T b ]\A®A> C N ~ k 

=■ C k ~ 2 MiC N ~ k (45) 

The operator Mj is calculated by expanding the commutator and using (^) and 
< ss\T B =< ss\, with the result 

Mi = < ss\l <8> f^A® A > - < ss\l ® f^A® A > 

= CM -CM (46) 

Using the fact that C = C — K,Ai = — X i; we rewrite Mj as CXi — KAi. Using the same 
argument as we did in the RS case we find 

<W\X lC ^\V> 

- < w\c«\v > 1 ' 

Similar manipulation shows that in the FS update the current is given by the same general 
form as in fl47|) , with C replaced with C. Writing C = C + K, expanding [C + KY N ~^ and 
using K\V >= 0, we find that the currents of each conserved species in these two updates 
are equal for arbitrary transition probabilities, i.e. 

<w\x i c N - 1 \v > , . 

< • r - >=< 3 - >= < w\c\v > (48 > 



We have proved this relation under a more general condition as in |19[ , thats we do not 
assmue C = C or equivalently K = 0. As we have pointed out in the introduction and will 
prove in the next subsection the condition K = puts physical restrictions on the steady 
state, namely in this case all the current correlators become distance- indepdent. Moreover 
we find the following relation between density profiles and currents in the two updates, which 
is valid regardless of the bulk and boundary transition rates: 

< ri >^ (k)- < n l >^ (k) =< >=< JV > (49) 
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the proof of which is easy to see, once we note that 



<n l >^{k) = J- < Wl&^A^-^V > 
Z N 

= -J- < W\C h -\Ai -Xi)C N ~ k \V > 
Z N 



< n l (k) - J% (50) 



This is the generalization of a similar relation for for 1-ASEP which is now valid for 
each species seperately and for arbitrary transition probabilities. 

3.1 Equal time current correlators 

In this section we consider only the RS update and find the MPA expressions for the current 
correlators < J % k J{ > of two kinds of particles i and j at sites k and 1. For I > k + 1 (i.e. 
disjoint links) we obtain a simple relation. Starting from the definition 

< JlJi >:=< S\[ht 1>k Mhf- ltl ,n]\P(t) > (51) 
and proceeding exactly along the lines which led to ([35"D we find 

< J k J i >n- <W \c»\v> (52) 

The proof of this relation is detailed in the appendix. For I = k + 1 (i.e. consecutive links) 

no simple relation is obtained. Several remarks are in order now. 

Remarks: 

1: In the special case K = 0, the two point correlator (and in fact all the n-point current 
correlators (see the appendix)), become site- independent, since in this case the operators X^ 
commute with C and the two point correlator can be written as 

< J k Jl >N- <W \ C N\y > 
2: For number- valued Xi and Xj this relation implies 

< J^jj > N = < J{J\ > N = < J* >N< J j >N-1 ■ (54) 

with similar relations for higher correlators. This is of course a finite-size effect and in the 
thermodynamic limit, the currents at disjoint links are not correlated. However this is true 
only when some of the X^s are number- valued . In this way we have shown that the nature 
of X, controlls the current correlators in a very definite way, namely for general X, but with 
K := YjXi = 0, the correlators are site independent and for number valued X, there is no 
correlation in the thermodynamic limit. 

3: No such simple relations can be obtained for the FS and BS updates even for disjoint 
links. This is due to that fact in these updates, no matter how far are the links, the cur- 
rent operators for these links contain a common string of local operators (see fl42|)). Thus 
although the currents of the two ordered updates are equal, their correlators may not be 
related to each other in any simple way. In this sense the steady state of the two updates 
may not be physically equivalent. 
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4 Beyond nearest neighbour interactions 

The matrix product ansatz as formulated in H, can be generalized to models with more 
general Hamiltonians.0. In the following we consider a Hamiltonian with nearest and next- 
nearest neighbor interactions, although our analysis can be generalized to more non-local 
Hamiltonians. Consider a Hamiltonian of the form: 

AT-2 

H = h 12 + h k,k+2 + h N -i,N- (55) 

k=l 

where h kk+2 acts on the three sites k, k + 1 and k + 2 and h\ 2 and h N _ 1N are boundary 
terms. Writing the steady state as in (||) one finds the following conditions for stationarity: 

h B A® A® A = X ® A- A® X (56) 

(h N . 1:N A®A-X)\V>=0 (57) 
< W\(h 12 A®A + X) = (58) 
where X is in general an operator valued tensor in h ® h, i.e. 

X :=Y J X a ^\a,(3> (59) 

and | a > and \(3 > stand for the states of one site. Denoting as before C :=< s\A > and 
K :=< ss\X > one finds again that equation (^) is true also in this case. To find the MPA 
expressions for the currents we proceed as in section 4 and find: 

j t < 4 >=< S\[h B k _ 2 , k + h B _ 1M1 + h B k+2 ,rl}\P(t) > (60) 
Local conservation of particles, now implies 

^kl = — [^f-l,fc+l> + ^fc+l] (61) 



Acting on by < sss\i l k <g)l®l + l(g)f£®l + l®l<g>f£ and using (|21|J22J) we find that 

[XF>+X®,C\ = 0, (62) 

where 



X^:=£Xy .Vr : (63) 



'(1) ._ V- y Y W ._ 

3 3 

Inserting ([H]) in (|60"D, we find again a continuity equation as in (B^) with the currents 

S -ii \i,B 



< 4 >=< S\[h B rl\ - [ht^rUWPit) > (64) 



3 The possibility of extending MPA to include more non-local interactions has also been recently noted in 
20| 
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The MPA expression for this current is obtained along the same line as in section 4. After 
similar manipulations we find 

J l k = -^C k -* <sss\[h B ,l®l®T i ]\A®A®A>C N - k 
Z N 

- -^C k - 2 < sssW^.r ®l®l]\A®A®A> C N - k+1 (65) 

Expanding the commutators and using ([H]) , (|2^) and (f)6"l) we find 

Ji = C k - 2 < ssKIBt' + t 4 ® 1)\X > C N - k (66) 

which finally yields 

p ^ <_ 

< W\C N \V > 



< j' >= < > (67 ) 



where 

A« + X„ (i ) (68) 

V 

The generalization of these results to Hamiltonians with more non-local interactions is ob- 
vious. 

5 Discussion 

We have discussed the general structure of MPA states in stochastic systems and have proved 
certain general relations for general stochastic processes in random and ordered updates. We 
have tried to be as general as possible, our results in sections 3-5 are independent of the 
bulk and boundary transition rates and are also independent of the asymmetry caused by 
driving. We have found general MPA expressions for the currents and current correlators 
and have shown that when in the MPA formalism one uses general operator-valued X but 
with J2i Xi = 0, the current correlators while non- vanishing, become site independent. 
And when any of the X^s are c- numbers, then the connected correlation functions of the 
corresponding currents vanish in the thermodynamic limit. 
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7 Appendix 

In this appendix we present the proof of (65). In the same spirit that we have derived the 
expressions for the currents we write 

<JlJi>=<W\X k Ji\V> (69) 
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where 



J k J l 



1 C k -M : C' h -M,C X 

N 

.( •••--!( ' \ _ /, A.\f "-'•--/ V , 

3 Ji 3 J 



-C k - 2 (CXi- KAi)C l - k - 2 (XjC - AjK)C N - 1 (70) 
Zn 



where in the last line we have used two equivalent expressions for M, in order to move the 
two K's to the left and right respectively and act by them on < W\ and \V > to obtain 

< J ^ >- < w|c^> 171 j 

This result is true for the two point correlators, since in higher correlators one can not 
eliminate all the K's in (|70f) . When = 0, this further simplifies to fl53f ). Furthermore in 
this case formula (p3|) trivialy generalizes to n-point functions. 
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